Hexahelicene DNA-binding: Minor groove selectivity, semi-intercalation and chiral recognition.

In this contribution, we focused on a fundamental study targeting the interaction of water-soluble [6]helicene derivative 1 (1-butyl-3-(2-methyl[6]helicenyl)-imidazolium bromide) with double-stranded (ds) DNA. A synthetic 30-base pair duplex, plasmid, chromosomal calf thymus and salmon DNA were investigated using electrochemistry, electrophoresis and spectroscopic tools supported by molecular dynamics (MD) and quantum mechanical approaches. Both experimental and theoretical work revealed the minor groove binding of 1 to the dsDNA. Both the positively charged imidazole ring and hydrophobic part of the side chain contributed to the accommodation of 1 into the dsDNA structure. Neither intercalation into the duplex DNA nor the stable binding of 1 to single-stranded DNA were found in topoisomerase relaxation experiments with structural components of 1, i.e. [6]helicene (2) and 1-butyl-3-methylimidazolium bromide (3), nor by theoretical calculations. Finally, the binding of optically pure enantiomers (P)-1 and (M)-1 was studied using circular dichroism spectroscopy, isothermal titration calorimetry and UV Resonance Raman (UVRR) methods. Using MD and quantum mechanical methods, minor groove and semi-intercalation were proposed for compound 1 as the predominant binding modes. From the UVRR findings, we also can conclude that 1 tends to preferentially interact with adenine and guanine residues in the structure of dsDNA.

Unique peptidic agonists of a juvenile hormone receptor with species-specific effects on insect development and reproduction.

Juvenile hormones (JHs) control insect metamorphosis and reproduction. JHs act through a receptor complex consisting of methoprene-tolerant (Met) and taiman (Tai) proteins to induce transcription of specific genes. Among chemically diverse synthetic JH mimics (juvenoids), some of which serve as insecticides, unique peptidic juvenoids stand out as being highly potent yet exquisitely selective to a specific family of true bugs. Their mode of action is unknown. Here we demonstrate that, like established JH receptor agonists, peptidic juvenoids act upon the JHR Met to halt metamorphosis in larvae of the linden bug, Pyrrhocoris apterus. Peptidic juvenoids induced ligand-dependent dimerization between Met and Tai proteins from P. apterus but, consistent with their selectivity, not from other insects. A cell-based split-luciferase system revealed that the Met-Tai complex assembled within minutes of agonist presence. To explore the potential of juvenoid peptides, we synthesized 120 new derivatives and tested them in Met-Tai interaction assays. While many substituents led to loss of activity, improved derivatives active at sub-nanomolar range outperformed hitherto existing peptidic and classical juvenoids including fenoxycarb. Their potency in inducing Met-Tai interaction corresponded with the capacity to block metamorphosis in P. apterus larvae and to stimulate oogenesis in reproductively arrested adult females. Molecular modeling demonstrated that the high potency correlates with high affinity. This is a result of malleability of the ligand-binding pocket of P. apterus Met that allows larger peptidic ligands to maximize their contact surface. Our data establish peptidic juvenoids as highly potent and species-selective novel JHR agonists.

Explaining the interaction of mangiferin with MMP-9 and NF-ƙß: a computational study.

Mangiferin is a glycosylated xanthone widely distributed in nature, which exhibits wide pharmacological activities, highlighting its anti-cancer properties. Mangiferin interferes with inflammation, lipid, and calcium signaling, which selectively inhibits multiple NFkB target genes as interleukin-6, tumor necrosis factor, plasminogen, and matrix metalloproteinase, among others. In this work, the interactions of this polyphenol with MMP-9 and NF-κß are characterized by using computational chemistry methods. The results show MMP-9 inhibition by mangiferina is characterized for the interact with the catalytic Zn atom through a penta-coordinate structure. It is also demonstrated through a strong charge transfer established between mangiferin and Zn in the QM/MM study. Concerning the mangiferin/NF-κß system, the 92.3% of interactions between p50 sub-unity and DNA are maintained with a binding energy of - 8.04 kcal/mol. These findings indicate that mangiferin blocks the p50-p65/DNA interaction resulting in the loss of the functions of this hetero-dimeric member and suggesting inhibition of the cancer progression. Experimental results concerning the anti-cancer properties of mangiferin show that this natural compound can inhibit selectively MMP-9 and NF-ƙß. Although the anti-tumor properties of mangiferin are well defined, its molecular mechanisms of actions are not described. In this work, a computational study is carried out to characterize the interactions of mangiferin with these molecular targets. The results obtained corroborate the anti-proliferative and anti-apoptotic activity of mangiferin and provide a depiction of its mechanisms of action.

Mechanistic insight into the RNA-stimulated ATPase activity of tick-borne encephalitis virus helicase.

The helicase domain of nonstructural protein 3 (NS3H) unwinds the double-stranded RNA replication intermediate in an ATP-dependent manner during the flavivirus life cycle. While the ATP hydrolysis mechanism of Dengue and Zika viruses NS3H has been extensively studied, little is known in the case of the tick-borne encephalitis virus NS3H. We demonstrate that ssRNA binds with nanomolar affinity to NS3H and strongly stimulates the ATP hydrolysis cycle, whereas ssDNA binds only weakly and inhibits ATPase activity in a noncompetitive manner. Thus, NS3H is an RNA-specific helicase, whereas DNA might act as an allosteric inhibitor. Using modeling, we explored plausible allosteric mechanisms by which ssDNA inhibits the ATPase via nonspecific binding in the vicinity of the active site and ATP repositioning. We captured several structural snapshots of key ATP hydrolysis stages using X-ray crystallography. One intermediate, in which the inorganic phosphate and ADP remained trapped inside the ATPase site after hydrolysis, suggests that inorganic phosphate release is the rate-limiting step. Using structure-guided modeling and molecular dynamics simulation, we identified putative RNA-binding residues and observed that the opening and closing of the ATP-binding site modulates RNA affinity. Site-directed mutagenesis of the conserved RNA-binding residues revealed that the allosteric activation of ATPase activity is primarily communicated via an arginine residue in domain 1. In summary, we characterized conformational changes associated with modulating RNA affinity and mapped allosteric communication between RNA-binding groove and ATPase site of tick-borne encephalitis virus helicase.

Specific anion effect on properties of HRV 3C protease.

Specific salts effect is intensively studied from the prospective of modification of different physico-chemical properties of biomacromolecules. Limited knowledge of the specific salts effect on enzymes led us to address the influence of five sodium anions: sulfate, phosphate, chloride, bromide, and perchlorate, on catalytic and conformational properties of human rhinovirus-14 (HRV) 3C protease. The enzyme conformation was monitored by circular dichroism spectrum (CD) and by tyrosines fluorescence. Stability and flexibility of the enzyme have been analyzed by CD in the far-UV region, differential scanning calorimetry and molecular dynamics simulations, respectively. We showed significant influence of the anions on the enzyme properties in accordance with the Hofmeister effect. The HRV 3C protease in the presence of kosmotropic anions, in contrast with chaotropic anions, exhibits increased stability, rigidity. Correlations of stabilization effect of anions on the enzyme with their charge density and the rate constant of the enzyme with the viscosity B-coefficients of anions suggest direct interaction of the anions with HRV 3C protease. The role of stabilization and decreased fluctuation of the polypeptide chain of HRV 3C protease on its activation in the presence of kosmotropic anions is discussed within the frame of the macromolecular rate theory.

Metallic Effects on p-Hydroxyphenyl Porphyrin Thin-Film-Based Planar Optical Waveguide Gas Sensor: Experimental and Computational Studies.

Metal effects on the gas sensing behavior of metal complexes of 5,10,15,20-tetrakis(4-hydroxyphenyl)porphyrin (THPP) thin film was investigated in terms of detecting NO2 gas by the planar optical waveguide. For this purpose, several THPP and metal complexes were synthesized with different central metal ions: Co(II), Ni(II), Cu(II), and Zn(II). Planar optical gas sensors were fabricated with the metalloporphyrins deposited on K+ ion-exchanged soda-lime glass substrate with the spin coating method serving as host matrices for gas interaction. All of the THPP complex's films were fully characterized by UV-Vis, IR and XPS spectroscopy, and the laser light source wavelength was selected at 520 and 670 nm. The results of the planar optical waveguide sensor show that the Zn-THPP complex exhibits the strongest response with the lowest detectable gas concentration of NO2 gas for both 520 nm and 670 nm. The Ni-THPP and Co-THPP complexes display good efficiency in the detection of NO2, while, on the other hand, Cu-THPP shows a very low interaction with NO2 gas, with only 50 ppm and 200 ppm detectable gas concentration for 520 nm and 670 nm, respectively. In addition, molecular dynamic simulations and quantum mechanical calculations were performed, proving to be coherent with the experimental results.

Interactions between a dsDNA Oligonucleotide and Imidazolium Chloride Ionic Liquids: Effect of Alkyl Chain Length, Part I.

Ionic liquids (ILs) have become nearly ubiquitous solvents and their interactions with biomolecules has been a focus of study. Here, we used the fluorescence emission of DAPI, a groove binding fluorophore, coupled with molecular dynamics (MD) simulations to report on interactions between imidazolium chloride ([Imn,1]+) ionic liquids and a synthetic DNA oligonucleotide composed entirely of T/A bases (7(TA)) to elucidate the effects ILs on a model DNA duplex. Spectral shifts on the order of 500-1000 cm-1, spectral broadening (~1000 cm-1), and excitation and emission intensity ratio changes combine to give evidence of an increased DAPI environment heterogeneity on added IL. Fluorescence lifetimes for DAPI/IL solutions yielded two time constants 0.15 ns (~80% to 60% contribution) and 2.36-2.71 ns for IL up to 250 mM. With DNA, three time constants were required that varied with added IL (0.33-0.15 ns (1-58% contribution), ~1.7-1.0 ns (~5% contribution), and 3.8-3.6 ns (94-39% contribution)). MD radial distribution functions revealed that π-π stacking interactions between the imidazolium ring were dominant at lower IL concentration and that electrostatic and hydrophobic interactions become more prominent as IL concentration increased. Alkyl chain alignment with DNA and IL-IL interactions also varied with IL. Collectively, our data showed that, at low IL concentration, IL was primarily bound to the DNA minor groove and with increased IL concentration the phosphate regions and major groove binding sites were also important contributors to the complete set of IL-DNA duplex interactions.

Surface Interaction of Ionic Liquids: Stabilization of Polyethylene Terephthalate-Degrading Enzymes in Solution.

The effect of aqueous solutions of selected ionic liquids solutions on Ideonella sakaiensis PETase with bis(2-hydroxyethyl) terephthalate (BHET) substrate were studied by means of molecular dynamics simulations in order to identify the possible effect of ionic liquids on the structure and dynamics of enzymatic Polyethylene terephthalate (PET) hydrolysis. The use of specific ionic liquids can potentially enhance the enzymatic hydrolyses of PET where these ionic liquids are known to partially dissolve PET. The aqueous solution of cholinium phosphate were found to have the smallest effect of the structure of PETase, and its interaction with (BHET) as substrate was comparable to that with the pure water. Thus, the cholinium phosphate was identified as possible candidate as ionic liquid co-solvent to study the enzymatic hydrolyses of PET.

Substituent Effect on Porphyrin Film-Gas Interaction by Optical Waveguide: Spectrum Analysis and Molecular Dynamic Simulation.

Substituent effect on optical gas sensing performance in porphyrin-based optical waveguide detection system was studied by molecular dynamics simulation (MDS), absorption/emission spectrum analysis, and optical waveguide (OWG) detection. The affinities of porphyrin with seven types of substituents (-H, -OH, -tBu, -COOH, -NH2, -OCH3, -SO3-) on para position of meso-phenyl porphyrin toward gas molecules in adsorption process were studied in different size of boxes with the same pressure and concentration. Analyte gases (CO2, H2S, HCl, NO2) were exposed to porphyrin film in absorption spectrophotometer, and in OWG with evanescent field excited by a guiding laser light with 670 nm wavelength. The extent of interaction between host molecule and the guest analytes was analyzed by the number of gas molecules in vicinity of 0.3 nm around substituents of porphyrin molecules. Optical waveguide results reveal that sulfonate porphyrin is mostly responsive to hydrochloride, hydrosulfide gas and nitrogen dioxide gases with strong response intensity. Molecular dynamics and spectral analysis provide objective information about the molecular state and sensing properties. Molecular rearrangements induced by gas exposure was studied by spectral analysis and surface morphology before and after gas exposure taking hydrosulfide gas as an example. Film-gas interaction mechanism was discussed in terms of each gas and substituent group characters.

Conserved Dynamic Mechanism of Allosteric Response to L-arg in Divergent Bacterial Arginine Repressors.

Hexameric arginine repressor, ArgR, is the feedback regulator of bacterial L-arginine regulons, and sensor of L-arg that controls transcription of genes for its synthesis and catabolism. Although ArgR function, as well as its secondary, tertiary, and quaternary structures, is essentially the same in E. coli and B. subtilis, the two proteins differ significantly in sequence, including residues implicated in the response to L-arg. Molecular dynamics simulations are used here to evaluate the behavior of intact B. subtilis ArgR with and without L-arg, and are compared with prior MD results for a domain fragment of E. coli ArgR. Relative to its crystal structure, B. subtilis ArgR in absence of L-arg undergoes a large-scale rotational shift of its trimeric subassemblies that is very similar to that observed in the E. coli protein, but the residues driving rotation have distinct secondary and tertiary structural locations, and a key residue that drives rotation in E. coli is missing in B. subtilis. The similarity of trimer rotation despite different driving residues suggests that a rotational shift between trimers is integral to ArgR function. This conclusion is supported by phylogenetic analysis of distant ArgR homologs reported here that indicates at least three major groups characterized by distinct sequence motifs but predicted to undergo a common rotational transition. The dynamic consequences of L-arg binding for transcriptional activation of intact ArgR are evaluated here for the first time in two-microsecond simulations of B. subtilis ArgR. L-arg binding to intact B. subtilis ArgR causes a significant further shift in the angle of rotation between trimers that causes the N-terminal DNA-binding domains lose their interactions with the C-terminal domains, and is likely the first step toward adopting DNA-binding-competent conformations. The results aid interpretation of crystal structures of ArgR and ArgR-DNA complexes.

The study of conformational changes in photosystem II during a charge separation.

Photosystem II (PSII) is a multi-subunit pigment-protein complex and is one of several protein assemblies that function cooperatively in photosynthesis in plants and cyanobacteria. As more structural data on PSII become available, new questions arise concerning the nature of the charge separation in PSII reaction center (RC). The crystal structure of PSII RC from cyanobacteria Thermosynechococcus vulcanus was selected for the computational study of conformational changes in photosystem II associated to the charge separation process. The parameterization of cofactors and lipids for classical MD simulation with Amber force field was performed. The parametrized complex of PSII was embedded in the lipid membrane for MD simulation with Amber in Gromacs. The conformational behavior of protein and the cofactors directly involved in the charge separation were studied by MD simulations and QM/MM calculations. This study identified the most likely mechanism of the proton-coupled reduction of plastoquinone QB. After the charge separation and the first electron transfer to QB, the system undergoes conformational change allowing the first proton transfer to QB- mediated via Ser264. After the second electron transfer to QBH, the system again adopts conformation allowing the second proton transfer to QBH-. The reduced QBH2 would then leave the binding pocket.

A residue of motif III positions the helicase domains of motor subunit HsdR in restriction-modification enzyme EcoR124I.

Type I restriction-modification enzymes differ significantly from the type II enzymes commonly used as molecular biology reagents. On hemi-methylated DNAs type I enzymes like the EcoR124I restriction-modification complex act as conventional adenine methylases at their specific target sequences, but unmethylated targets induce them to translocate thousands of base pairs through the stationary enzyme before cleaving distant sites nonspecifically. EcoR124I is a superfamily 2 DEAD-box helicase like eukaryotic double-strand DNA translocase Rad54, with two RecA-like helicase domains and seven characteristic sequence motifs that are implicated in translocation. In Rad54 a so-called extended region adjacent to motif III is involved in ATPase activity. Although the EcoR124I extended region bears sequence and structural similarities with Rad54, it does not influence ATPase or restriction activity as shown in this work, but mutagenesis of the conserved glycine residue of its motif III does alter ATPase and DNA cleavage activity. Through the lens of molecular dynamics, a full model of HsdR of EcoR124I based on available crystal structures allowed interpretation of functional effects of mutants in motif III and its extended region. The results indicate that the conserved glycine residue of motif III has a role in positioning the two helicase domains.

Aggregation and metal-complexation behaviour of THPP porphyrin in ethanol/water solutions as function of pH.

The effect of pH change on 5,10,15,20-Tetrakis(4-hydroxyphenyl)-21H,23H-porphine (THPP) with its aggregation as function of water-ethanol mixture was studied with UV-vis, fluorescence, Raman and computational analysis. In neutral pH, THPP was present as free-base and, increasing the water amount, aggregation occurred with the formation of H- and J-aggregates. The aggregation constant and the concentration of dimers were calculated, other information about the dimer aggregation were evaluated by computational study. In acidic pH, by the insertions of two hydrogens in the porphyrin rings, the porphyrin changed its geometry with a ring deformation confirmed by red-shifted spectrum and quenching in fluorescence; at this low pH, increasing the water amount, the acidic form (THPPH2)2+ resulted more stable due to a polar environment with stronger interaction by hydrogen bonding. In basic pH, reached by NH4OH, THPP porphyrin was able to react with alkali metals in order to form sitting-atop complex (M2THPP) confirmed by the typical absorption spectrum of metallo-porphyrin, Raman spectroscopy and by computational analysis.

Transmembrane helix connectivity in Orai1 controls two gates for calcium-dependent transcription.

The channel Orai1 requires Ca2+ store depletion in the endoplasmic reticulum and an interaction with the Ca2+ sensor STIM1 to mediate Ca2+ signaling. Alterations in Orai1-mediated Ca2+ influx have been linked to several pathological conditions including immunodeficiency, tubular myopathy, and cancer. We screened large-scale cancer genomics data sets for dysfunctional Orai1 mutants. Five of the identified Orai1 mutations resulted in constitutively active gating and transcriptional activation. Our analysis showed that certain Orai1 mutations were clustered in the transmembrane 2 helix surrounding the pore, which is a trigger site for Orai1 channel gating. Analysis of the constitutively open Orai1 mutant channels revealed two fundamental gates that enabled Ca2+ influx: Arginine side chains were displaced so they no longer blocked the pore, and a chain of water molecules formed in the hydrophobic pore region. Together, these results enabled us to identify a cluster of Orai1 mutations that trigger Ca2+ permeation associated with gene transcription and provide a gating mechanism for Orai1.

Theoretical and experimental study of the antifreeze protein AFP752, trehalose and dimethyl sulfoxide cryoprotection mechanism: correlation with cryopreserved cell viability.

In this work the physico-chemical properties of selected cryoprotectants (antifreeze protein TrxA-AFP752, trehalose and dimethyl sulfoxide) were correlated with their impact on the constitution of ice and influence on frozen/thawed cell viability. The freezing processes and states of investigated materials solutions were described and explained from a fundamental point of view using ab-initio modelling (molecular dynamics, DFT), Raman spectroscopy, Differential Scanning Calorimetry and X-Ray Diffraction. For the first time, in this work we correlated the microscopic view (modelling) with the description of the frozen solution states and put these results in the context of human skin fibroblast viability after freezing and thawing. DMSO and AFP had different impacts on their solution's freezing process but in both cases the ice crystallinity size was considerably reduced. DMSO and AFP treatment in different ways improved the viability of frozen/thawed cells.

The mechanism of the glycosylase reaction with hOGG1 base-excision repair enzyme: concerted effect of Lys249 and Asp268 during excision of 8-oxoguanine.

The excision of 8-oxoguanine (oxoG) by the human 8-oxoguanine DNA glycosylase 1 (hOGG1) base-excision repair enzyme was studied by using the QM/MM (M06-2X/6-31G(d,p):OPLS2005) calculation method and nuclear magnetic resonance (NMR) spectroscopy. The calculated glycosylase reaction included excision of the oxoG base, formation of Lys249-ribose enzyme-substrate covalent adduct and formation of a Schiff base. The formation of a Schiff base with ΔG# = 17.7 kcal/mol was the rate-limiting step of the reaction. The excision of the oxoG base with ΔG# = 16.1 kcal/mol proceeded via substitution of the C1΄-N9 N-glycosidic bond with an H-N9 bond where the negative charge on the oxoG base and the positive charge on the ribose were compensated in a concerted manner by NH3+(Lys249) and CO2-(Asp268), respectively. The effect of Asp268 on the oxoG excision was demonstrated with 1H NMR for WT hOGG1 and the hOGG1(D268N) mutant: the excision of oxoG was notably suppressed when Asp268 was mutated to Asn. The loss of the base-excision function was rationalized with QM/MM calculations and Asp268 was confirmed as the electrostatic stabilizer of ribose oxocarbenium through the initial base-excision step of DNA repair. The NMR experiments and QM/MM calculations consistently illustrated the base-excision reaction operated by hOGG1.

The helical domain of the EcoR124I motor subunit participates in ATPase activity and dsDNA translocation.

Type I restriction-modification enzymes are multisubunit, multifunctional molecular machines that recognize specific DNA target sequences, and their multisubunit organization underlies their multifunctionality. EcoR124I is the archetype of Type I restriction-modification family IC and is composed of three subunit types: HsdS, HsdM, and HsdR. DNA cleavage and ATP-dependent DNA translocation activities are housed in the distinct domains of the endonuclease/motor subunit HsdR. Because the multiple functions are integrated in this large subunit of 1,038 residues, a large number of interdomain contacts might be expected. The crystal structure of EcoR124I HsdR reveals a surprisingly sparse number of contacts between helicase domain 2 and the C-terminal helical domain that is thought to be involved in assembly with HsdM. Only two potential hydrogen-bonding contacts are found in a very small contact region. In the present work, the relevance of these two potential hydrogen-bonding interactions for the multiple activities of EcoR124I is evaluated by analysing mutant enzymes using in vivo and in vitro experiments. Molecular dynamics simulations are employed to provide structural interpretation of the functional data. The results indicate that the helical C-terminal domain is involved in the DNA translocation, cleavage, and ATPase activities of HsdR, and a role in controlling those activities is suggested.

Quantum Calculations Indicate Effective Electron Transfer between FMN and Benzoquinone in a New Crystal Structure of Escherichia coli WrbA.

UNLABELLED: Quantum mechanical calculations using the Marcus equation are applied to compare the electron-transfer probability for two distinct crystal structures of the Escherichia coli protein WrbA, an FMN-dependent NAD(P)H: quinone oxidoreductase, with the bound substrate benzoquinone. The calculations indicate that the position of benzoquinone in a new structure reported here and solved at 1.33 Å resolution is more likely to be relevant for the physiological reaction of WrbA than a previously reported crystal structure in which benzoquinone is shifted by ∼5 Å. Because the true electron-acceptor substrate for WrbA is not yet known, the present results can serve to constrain computational docking attempts with potential substrates that may aid in identifying the natural substrate(s) and physiological role(s) of this enzyme. The approach used here highlights a role for quantum mechanical calculations in the interpretation of protein crystal structures.

Influence of ligand binding on structure and thermostability of human α1-acid glycoprotein.

Ligand binding of neutral progesterone, basic propranolol, and acidic warfarin to human α1-acid glycoprotein (AGP) was investigated by Raman spectroscopy. The binding itself is characterized by a uniform conformational shift in which a tryptophan residue is involved. Slight differences corresponding to different contacts of the individual ligands inside the ß-barrel are described. Results are compared with in silico ligand docking into the available crystal structure of deglycosylated AGP using quantum/molecular mechanics. Calculated binding energies are -18.2, -14.5, and -11.5 kcal/mol for warfarin, propranolol, and progesterone, respectively. These calculations are consistent with Raman difference spectroscopy; nevertheless, minor discrepancies in the precise positions of the ligands point to structural differences between deglycosylated and native AGP. Thermal dynamics of AGP with/without bounded warfarin was followed by Raman spectroscopy in a temperature range of 10-95 °C and analyzed by principal component analysis. With increasing temperature, a slight decrease of α-helical content is observed that coincides with an increase in ß-sheet content. Above 45 °C, also ß-strands tend to unfold, and the observed decrease in ß-sheet coincides with an increase of ß-turns accompanied by a conformational shift of the nearby disulfide bridge from high-energy trans-gauche-trans to more relaxed gauche-gauche-trans. This major rearrangement in the vicinity of the bridge is not only characterized by unfolding of the ß-sheet but also by subsequent ligand release. Hereby, ligand binding alters the protein dynamics, and the more rigid protein-ligand complex shows an improved thermal stability, a finding that contributes to the reported chaperone-like function of AGP.

Molecular dynamics comparison of E. coli WrbA apoprotein and holoprotein.

WrbA is a novel multimeric flavodoxin-like protein of unknown function. A recent high-resolution X-ray crystal structure of E. coli WrbA holoprotein revealed a methionine sulfoxide residue with full occupancy in the FMN-binding site, a finding that was confirmed by mass spectrometry. In an effort to evaluate whether methionine sulfoxide may have a role in WrbA function, the present analyses were undertaken using molecular dynamics simulations in combination with further mass spectrometry of the protein. Methionine sulfoxide formation upon reconstitution of purified apoWrbA with oxidized FMN is fast as judged by kinetic mass spectrometry, being complete in ∼5 h and resulting in complete conversion at the active-site methionine with minor extents of conversion at heterogeneous second sites. Analysis of methionine oxidation states during purification of holoWrbA from bacterial cells reveals that methionine is not oxidized prior to reconstitution, indicating that methionine sulfoxide is unlikely to be relevant to the function of WrbA in vivo. Although the simulation results, the first reported for WrbA, led to no hypotheses about the role of methionine sulfoxide that could be tested experimentally, they elucidated the origins of the two major differences between apo- and holoWrbA crystal structures, an alteration of inter-subunit distance and a rotational shift within the tetrameric assembly.